r] transitions between charmonia with meson loop contributions 
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We study the r\ transitions between i/r(4040/4160) and J /if/ by introducing charmed meson loops in an effective 
Lagrangian approach to enhance the decay amplitudes. The branching fractions S[i/f(4040) — » J/if/rfi and 
Sh//(4160) — > J/if/rj] estimated in this paper can remarkably explain the experimental measurements of Belle 
and BESIII within a reasonable parameter range. The rf transition between i/r(4160) and J/if/ is also investigated, 
and the branching fraction is under the upper limit of CLEO, which can be tested by future experiments. 
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I. INTRODUCTION 

Among many heavy quarknonia, the chamonium sector es- 
pecially has much abundant spectroscopy and decay modes 
observed [ 1 ] . In their mass range, there are also many discrep- 
ancies between theoretical predictions and experimental mea- 
surements because of the non-pertubative property of Quan- 
tum Chromodynamics (QCD). The study of spectroscopy 
and decay behavior of the charmonia undoubtedly enrich our 
knowledge of how QCD works in hadron physics. 

When considering the decay process ty\ — > if/iP, the di- 
rect coupling among the initial and final charmonia, if/\ and 
t/>2, and a light chiral meson, f, is highly suppressed due to 
the OZI rule. Because Belle and BESIII have recently ob- 
served branching fractions for some processes of this kind, of 
the order of 1CT 3 , hence we need to consider other mechanism 
which enhances this type of decay amplitude. The charmonia, 
whose masses are above the threshold of a charmed meson 
pair, dominantly decay into a pair of charmed mesons, which 
can couple with a light meson and a charmonium by exchang- 
ing a proper charmed meson in the final states. The contri- 
bution from such a structure, i.e., meson loop contribution, is 
dominant in this decay and becomes important in understand- 
ing decay behaviors of higher charmonia. 

Taking (^(3770) as an example, which is the first charmo- 
nium above the threshold of open charmed mesons, the BES 
Collaboration announced that the branching fraction of its 
non-DD decay is £[^(3770) -> non-DD] = (14.7±3.2)% [2- 
5]. Such a large non-DD branching fraction is several times 
larger than expected in theory [6]. To resolve this discrepancy, 
the authors in Refs. [7, 8] have taken account of the charmed 
meson loop, where the initial charmonium 013770) decays 
into a charmed meson pair DD and this pair couples to a 
vector and/or pseudoscalar meson by exchanging a D meson. 
After including the meson loop contributions, the large non- 
DD branching fraction of i/r(3770) can be nicely explained. 
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The decay mode i/r(3770) — > J/0 3 and the lineshape around 
0(3770) have been studied with meson loop contributions in 
Ref. [9]. 

Even though the measurements of the higher charmonia are 
not yet enough to discuss rj transitions, there appear some 
experiments of this kind. In the International Conference 
on High Energy Physics, the Belle Collaboration reported 
their measurements for 77 transitions between 0(4160/4040) 
and J/i// [10]. They announced that £[0(4040) -> rjj/if/] ■ 
T e + e - (0(4040)) = 5.1 ±0.8 ± 1.1 eV or 12.4 ± 1.2 ± 1.2 eV with 
different fitting parameters to the data. The corresponding re- 
sults for *A(4160) are £[0(4160) -» 7] J /iff] ■ r e+e - (0(4160)) = 
4. 1 ±0.5±0.8 eV or 15.2+1 .2+ 1 .5 eV. Taking T e+ir (0(4040)) = 
(0.86 + 0.07) keV and IW(0(4 160)) = (0.83 + 0.07) keV, one 
obtains the branching ratios as £[0(4040) — > J/if/rj] = (0.59 + 
0.11 ± 0.14)% or (1.44 + 0.18 + 0.18)% and £[0(4160) -> 
J/i/rrj] = (0.50 ± 0.07 + 0.11)% or (1.83 ± 0.21 ± 0.24)%. 
Before this measurement, only the upper limits of branching 
ratios for 0(4040) — » rfj/i// and 0(4160) — > rfj/i// were re- 
ported by the CLEO Collaboration, which are < 7 x 10~ 3 and 
< 8 x 10 3 [1 1], respectively. Recently, the BESIII Collabora- 
tion also analyzed the production of e + e~ — > r/J/if/ at a center- 
of-mass energy of yfs = 4.009 GeV. Because the Born cross 
section is reported to be (32.1 + 2.8 ± 1.3) pb [12], the corre- 
sponding fractional transition rate is £[i^(4040) — > rjJ/iff] = 
(5.2 + 0.5 ± 0.2 + 0.5) x 10~ 3 , which is consistent with the 
first solution of the Belle Collaboration and measurement by 
CLEO. 

Similar to the case of ^(3770), 1^(4040) and ^(4160) are 
above the threshold of charmed meson pairs, and dominantly 
decay into these. The experimental measurements stimulate 
us to study the 77 transition between 1^(4040/4160) and J /if/ 
with the meson loop mechanism, which is essential to under- 
stand the hidden charm decay behavior like higher charmonia. 

This paper is organized as follows. After introduction, a 
brief review of meson loop mechanism is presented and the 
corresponding amplitudes are calculated using an effective 
Lagrangian in Section II. Our numerical results of the branch- 
ing ratios are given in Section III. Section IV is devoted to 
summary. 
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H. 77 TRANSITION OF i//(4040) AND ifr(4 1 60) 

The meson loop effect plays a crucial role in understand- 
ing the 77 transition between the higher charmonium and // iff. 
Taking (^(4040) — > J/ifrrj as an example, involvement of the 
meson loop in the decay is depicted in Fig. 1. The charmo- 
nium i^(4040) is decomposed into £) ( * ) Z) ( * ) and by exchanging 
D or D* meson, the charmed meson pair converts itself into 
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FIG. 1 : The typical meson loop diagrams contributing to 1^(4040) — > 
J/iftr]. The initial <^(4040) couples with a charmed meson pair 
£)(*)£)W ; an( j 5y exchanging the D meson (left) or D* meson (right), 
the charmed meson pair converts itself into J/if/r] in the final state. 

The effective Lagrangian approach is adopted to evaluate 
the meson loop contributions to higher charmonia decay into 
J/if/T] as shown in Fig. 1. Utilizing the heavy quark limit and 
chiral symmetry, the involved effective Lagrangians read as 
[13, 14]: 
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with £> w = (D ( * )0 , D ( 5 * )+ ). Considering 77 and 77' mixing, 
one has P in the form, 
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and we adopt 9 = -19.1° in the present work [15, 16]. 

In the symmetric limit, the coupling constants between 
J/iff and charmed mesons satisfy g^DD = g^o'D'tiD/mo* = 
g^D'D y/m D m D - = mj/^/fj/y with to 7/iA and fj^ being the 
mass and decay constant of Jjif/. On the other hand, the 



coupling between pseudoscalar meson and charmed mesons 
can be related to the gauge coupling constant g by gp^'D" = 
g-PD*Dl ^m D m D > = 2g/f n with g = 0.59 and f n = 132 
MeV. The coupling constants between higher charmonia, such 
as i/r(4040) and 014 160), and charmed mesons are evaluated 
by the partial decay width assuming that these two char- 
monia dominantly decay into DD, D*D + h.c, and D*D* 
due to the restriction of the phase space. Furthermore 
the BaBar Collaboration has measured the ratios between 
these decay modes for i/r(4040) and i^(4160) [17], which are 
Br(i/<4040) -> DD)/Br(0(4O4O) -» D*D) = 0.24 + 0.05 ± 
0.12, Br( i^(4040) — > D*D*)/Br(i/f(4040) — > D*D) = 0.18 ± 
0.14 + 0.03 and Br(0(416O) -> DD)/Br(0(416O) -» D*D*) = 
0.02, Br(tA(4160) -> D*D)/Br(i/r(4160) -> D*D") = 0.34. 
Given the total decay width F^(4040) = 80 MeV and F^(4i60) = 
103 MeV, one obtains |§i/,(4040)ddI — 2.13, |^(4040)d*dI — 
1.70 GeV |^(4040)d*d*I = 3.34, |^ (4 i60)z)dI = 0.57, 
I^(4160)o«dI = 0.77 GeV" 1 , and \g^4m)D'D'\ = 2.23. 

With the Lagrangians listed above, we can obtain 
the hadronic decay amplitudes for i/r(4040)(;?o) >— » 
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where the amplitude M corresponds to the process in 
which the initial i/r(4040) is decomposed into a charmed me- 
son pair M\M 2 by exchanging M3, and then this meson pair 
converts itself into J/tfrrj in the final state. The form factor 
T{q 2 ,m 2 E ) = (m| - A 2 )/(q 2 - A 2 ) is introduced to evade 
infrared divergence in the loop integrals as well as to de- 
scribe the structure effects and off shell influence of the ex- 
changed mesons. The parameter A can be reparameterized 
as A = me + uKqcd with Aqcd = 220MeV, and the unique 
parameter a is taken to be of the order of unity [18]. 

To estimate the branching ratios of if/ — > J/ifrrj {if/' - 
{i/f(4040), (/'(4160)}) derived from hadronic loops, we summa- 
rize the amplitudes as, 
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where M\, M2 and AI3 are the amplitudes corresponding 
to the channels DD, D*D + h.c. and D*D*, respectively. 
The amplitudes can be reduced to a simple Lorentz struc- 
ture Sfwafi^, tj/^p"/^!^ multiplied with the coupling constants 



^ , which can be evaluated by the loop integrals. 
Since we can obtain only the absolute values of the coupling 
between (^(4040/4 160) and charmed mesons, we should, in 
general, introduce phase angles among different amplitudes. 
The total amplitudes can be expressed as 



Mot =Mi + e'* 1 A^2 + e'^-Mi, 



and the partial decay width is 

1 1 \Pn\ 



i JOT m^, 
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where the overline indicates the sum over the polarization vec- 
tors of iff' and J /if/ and \p n \ = A 1/2 (m 2 J ,,m 2 J/i// ,m 2 )/(2m^ l ) with 

A(x, y, z) - x 2 +y 2 + z 2 - 2xy - 2xz - zyz . 
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FIG. 2: (color online) The a dependence of the absolute values of 
coupling constants derived from meson loop contributions to 

1^(4040) — > J/if/ij (upper panel) and i^(4160) — > J/if/rj (lower panel). 
Different intermediate channels are considered; the green dotted, red 
dashed and black solid curves are corresponding to contributions 
from D'D*, D'D + h.c. and DD intermediates, respectively. 



III. NUMERICAL RESULTS AND DISCUSSION 

The coupling constant g®,*j21 defined in Eq. (5) can be 
evaluated by estimating the corresponding loop integrals. In 
general, the evaluated coupling constants from meson loop 
contributions are complex. In Fig. 2, we show the absolute 
values of the coupling constants corresponding to different in- 
termediate channels. For 1^(4040) — > J/if/rj process, the cou- 
pling constant corresponding to D*D + h.c channel is about 3 
times larger than that for DD channel, while the one for D*D* 
is one order larger than that for DD channel. The coupling 
constants for ifr(4160)J / ifrr/ are very similar to those for the 
1^(4040) — » //(/"7- The dominant contribution comes from the 
process with D*D* channel and the one corresponding to DD 
is very small and can be neglected. 

With the coupling evaluated above, one can estimate the 
branching ratios of 0-(4O4O/416O) — > J/ifnj. The phase angles 
in Eq. (6) are kept unknown, and the sign of the coupling 
constant for i//D w D ( *>, iff' = {^(4040), tA(4160)} can also be 
absorbed by these phase angles. Since the contributions from 
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FIG. 3: (color online) The a dependence of the branching ratios for 
0(4040) — > J/if/r] (upper panel) and 0(4160) — > J/if/r/ (lower panel) 
derived from meson loop. The black solid lines are the results corre- 
sponding to (f> = 0, and the red dashed lines are evaluated under the 
assumption = n. For comparison, the experimental measurements 
are also presented as cyan and yellow bands for different solutions 
to the data fitting [10] and blue band for measurements of BESIII 
Collaboration [12]. 
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FIG. 4: (color online) The a dependence of the coupling constants 
(upper panel) and the branching ratio (lower panel) of 0(4160) — » 
J/if/r]'. The coupling constants derived from different intermediate 
channels are presented and the branching ratios with different (f>2 = 
and <p2 = n correspond to the solid and dashed curves in the lower 
panel, respectively. 



DD channel can be negligible as indicated by the coupling 
constants shown in Fig. 2. For simplicity, we just fix the 
phase angle between M\ and M2 to be <p\ = 0. As for 02, in 
the present work we consider two extreme cases, 02 — and 
<p2 = n, which are corresponding to the black solid and red 
dashed curves in Fig. 3. 

For comparison, we also present the experimental measure- 
ments of Belle [10] marked as I and II for different solutions 
and results of BESIII [12] marked as BES with different hor- 
izontal bands in Fig. 3. Since the second solutions of the 
Belle Collaboration is not consistent with the one of the BE- 
SIII and CLEO Collaboration, we consider only the first so- 
lutions from the Belle Collaboration and results of BESIII 
Collaboration, which are the cyan and blue bands in Fig. 3. 
Our theoretical curves with $2= n have a reasonable overlap 
with the experimental measurements of the Belle Collabora- 
tion, which are about 0.7 ~ 1 .0 both for 0(4040) -> Jj\pj] and 



0(4160) — > J/if'T], while the overlap between theoretical curve 
and experimental data of the BESIII Collaboration is in the 
region a = 0.7 ~ 0.9. As for the case <p2 = 0, the overlaps be- 
tween our results and the experimental measurements appear 
in a relatively small a region. However, the values of a are 
still in a reasonable range. 

Other than the rj transition between 0(4040/4160) and J /if/, 
we also estimate the meson loop contributions to 0(4160) — * 
7/077'. The coupling constants derived from meson loops 
with different intermediate channels are presented in the up- 
per panel of Fig. 4. Similar to the 77 transitions between 
0(4040/4160) and J/ip, D*D* intermediate channel is dom- 
inant and the process via DD is negligible. Fixing (f>\ = 0, we 
can get two theoretical curves with (p2 - and 02 = n, which 
correspond to the black solid and red dashed lines in the lower 
panel of Fig. 4, respectively. In the region 0.7 < a < 0.9, the 
branching ratios are (4.3 ~ 6.8) x 10~ 4 and (2.2 ~ 3.6) x 10~ 4 
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corresponding to 4>2 - 0, n, respectively. These results are 
consistent with the CLEO measurement < 5 x 10~ 3 with 90% 
confidence level [11]. 

IV. SUMMARY 



The theoretical estimates have shown overlaps with the exper- 
imental measurements in a reasonable parameter range. More 
over, we have predicted the branching ratio of ^(4160) — > 
J/if/T]' of the order of 10 4 , which can be tested by future ex- 
periments. 



The charmonium above the threshold of a pair of charmed 
mesons dominantly decays into charmed mesons, which can 
couple with charmonium and light meson by exchanging a 
proper charmed meson, such a mechanism, i.e, the meson loop 
effect, is essential to understand the decay behavior of higher 
charmonia above the thresholds in obtaining the enhanced de- 
cay amplitudes. 

Stimulated by the experimental measurements by the Belle 
and BESIII Collaboration [10, 12], we introduce the me- 
son loop mechanism to study the rj transitions between 
^(4040/4160) and J/if/ in an effective Lagrangian approach. 
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